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Abstract Electron capture and positron decay rates 
are calculated for neutron-deficient Kr and Sr waiting 
point nuclei in stellar matter. The calculation is per- 
formed within the framework of pn-QRPA model for 
rp-process conditions. Fine tuning of particle-particle, 
particle-hole interaction parameters and a proper choice 
of the deformation parameter resulted in an accurate re- 
production of the measured half-lives. The same model 
parameters were used to calculate stellar rates. Inclu- 
sion of measured Gamow- Teller strength distributions 
finally led to a reliable calculation of weak rates that 
reproduced the measured half-lives well under limiting 
conditions. For the rp-process conditions, electron cap- 
ture and positron decay rates on ^^Kr and ^^Sr are of 
comparable magnitude whereas electron capture rates 
on ^^Sr and ^^Kr are 1-2 orders of magnitude bigger 
than the corresponding positron decay rates. The pn- 
QRPA calculated electron capture rates on ^''^Kr are 
bigger than previously calculated. The present calcu- 
lation strongly suggests that, under rp-process condi- 
tions, electron capture rates form an integral part of 
weak-interaction mediated rates and should not be ne- 
glected in nuclear reaction network calculations as done 
previously. 
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1 Introduction 

At low temperatures (Tg < 0.3K), energy generation 
and accompanying nucleosynthesis in explosive hydro- 
gen burning are characterized by the hot CNO cycles 
and at higher temperatures by the rp- and the ap- 
Tocesses (see for example (jll)). Recent studies (e.g. 

) have pointed toward extreme hydrogen burning 
where at sufficiently high temperature (Tg < 0.8K) and 
density conditions {p > 10'* gcm^'^), depending on the 
time scale of the explosive event, the rp- and the ap- 
processes reaction path may well proceed beyond mass 
A =64 and Z = 32. 

The rp-process is characterized by proton capture 
reaction rates that are orders of magnitude faster than 
any other competing process, specially /3-decay rates. 
The reaction path follows a series of fast (p,7)-reactions 
until further proton capture is inhibited, either by neg- 
ative proton capture Q-values (proton decay) or small 
positive proton capture Q-values (photodisintegration) . 
Now the reaction flow has to wait for the relatively 
slow /3-decay and the respective nucleus is hence called 
a waiting point. It is primarily the reliable estimate 
of half-lives of the waiting point nuclei that ultimately 
determine the timescale of the nucleosynthesis process 
and the produced isotopic abundances. The rp-process 
produces rapid nucleosynthesis on the proton-rich side 
of stability line toward heavier proton-rich nuclei. Such 
processes typically occur when hydrogen fuel is ignited 
under highly degenerate conditions in explosive events 
on the surface of compact objects like white dwarfs (no- 
vae) and neutron stars (X-ray bursts). 

Type I X-ray bursts have been pointed out as pos- 
sible sites for high temperature hydrogen burning via 
the rp-process (JSJ). The X-ray bursts are typically gen- 
erated by a thermonuclear runaway in the hydrogen- 
rich environment of an accreting compact object fed 
from a binary companion. The timescale for the ther- 
mal runaway ranges from 10-100 s ^. Within this 



timescale, the rp-process may proceed well beyond ^^Ni 
(|j). Previous attenipts to extend the network calcula- 
tions beyond Z = 32 suffered considerably due to lack 
of reliable nuclear physics data (JSJ). It was pointed out 
by authors in Rcf. ^ that important parameters for 
a successful rp-process nucleosynthesis calculation in- 
clude nuclear masses, nuclear deformations (specially 
in the regime A — 70-80 because of the wide variety 
of nuclear shapes displayed in the region) and finally 
the reliable calculation of stellar electron capture and 
/3-decay rates of the waiting point nuclei along a given 
reaction path since they determine time structure and 
abundance patterns. 

This article reports on the microscopic calculation 
of weak-interaction mediated rates (electron capture 
and /3-decay rates) on waiting point nuclei ■^2,74j^j. ^^j^j 
76,78gj._ fj^ jg |.p ^g noted that in a true sense only the 
even-even N ~ Z nuclei, ^^Kr and ^^Sr, act as waiting 
point nuclei because the next nucleus in the proton- 
capture chain is proton unbound. However, at times, 
the deformed isotopes, ^^Kr and ^^Sr, are also consid- 
ered as waiting point nuclei in literature (see e.g. (jj)). 
It is in this sense that 'i'^^'^'^Kr and 'i'6,78gj. g^^.^ iQYmed as 
waiting point nuclei in this manuscript.) The nuclear 
model chosen to perform this task is the proton-neutron 
quasiparticle random phase approximation (pn-QRPA) 
which has a proven track-record for the calculation of 
electron capture and /3-decay rates. Half-lives of f3^ 
decays were calculated systematically for about 6000 
neutron-rich nuclei between the beta stability line and 
the neutron drip line using the pn-QRPA model (l5|). 
Similarly half- lives for /3+/EC (electron capture) decays 
for neutron-deficient nuclei with atomic numbers Z = 
10 - 108 were calculated up to the proton drip line for 
more than 2000 nuclei using the same model (|6J). These 
microscopic calculations gave a remarkably good agree- 
ment with the then existing experimental data (within 
a factor of two for more than 90% (73%) of nuclei with 
experimental half-lives shorter than 1 s for /3~ (/3+/EC) 
decays). Most nuclei of interest of astrophysical impor- 
tance are the ones far from stability and one has to 
rely on theoretical models to estimate their beta de- 
cay properties. The accuracy of the pn-QRPA model 
increases with increasing distance from the /3-stability 
line (shorter half- lives) (l5|;|6|). This is a promising fea- 
ture with respect to the prediction of experimentally 
unknown half-lives (specially those present in the stel- 
lar interior), implying that the predictions are made 
on the basis of a realistic physical model. Later Nabi 
and Klapdor-Kleingrothaus reported the calculation of 
weak-interaction rates for more than 700 nuclei with 
A = 18 to 100 in stellar environment using the same 
nuclear model ([tI). These included capture rates, de- 
cay rates, gamma heating rates, neutrino energy loss 



rates, probabilities of beta-delayed particle emissions 
and energy rate of these particle emissions. The au- 
thors then presented a detailed calculation of stellar 
weak-interaction rates over a wide range of temperature 
and density scale for sd- dsj) and fp/fpg-shell nuclei (|9|) . 
Since then these calculations were further refined with 
use of more efficient algorithms, computing power, in- 
corporation of latest data from mass compilations and 
experimental values, and fine-tuning of model parame- 
ters both in the sd- shell (e.g. (jlOh) and fp-shell nuclei 

(e.g. ini)). 

The next section discusses the calculation of stel- 
lar weak-interaction mediated rates using the pn-QRPA 
formalism. The calculated Gamow- Teller strength dis- 
tributions of waiting point nuclei ^^'^''Kr and "^^'"^^Sr 
are also presented in this section. The rp-process elec- 
tron capture and positron decay rates are presented and 
compared with previous calculations in Section 3. Sum- 
mary and conclusions finally follow in Section 4. 



2 Description of the microscopic model 

The microscopic calculation of Gamow- Teller strength 
distributions under stellar conditions is a challenging 
task. The calculations become more challenging as the 
number of nucleons increases. Out of the very lim- 
ited available options, the pn-QRPA model can handle 
any arbitrarily heavy system of nucleons and performs 
a microscopic calculation of ground and excited state 
Gamow- Teller strength functions (|l2l:ll3l). 

For the present pn-QRPA calculation single parti- 
cle energies and wave functions were calculated in the 
Nilsson model (J14l ). which takes into account nuclear 
deformations. Pairing was treated in the BCS approx- 
imation. The proton-neutron residual interactions oc- 
cur in two different forms, namely as particle-hole and 
particle-particle interactions. These interactions were 
given separable form and were characterized by two in- 
teraction constants x f-nd k, respectively. These pa- 
rameters arc regarded as the two most important model 
parameters in the pn-QRPA theory (see Refs. (JSJ; l6J)). 
In this work a fine tuning of these parameters was done 
which optimally reproduced the measured half-lives for 
72,74j^j. and 'i'ejsgj, j-jTqj, measured half-lives sec Refs. 
( 15t lid : Il7|)). The particle- hole-interaction parameter, 
X, is known to infiuence the calculated position of the 
Gamow- Teller giant resonance (GTGR). Normally for 
the case of stable nuclei, the constant x is parame- 
terized such that the experimentally observed energy 
of the GTGR is reproduced in the best way. How- 
ever Hirsch and collaborators (|6|), using the pn-QRPA 
model, deduced the locally best x values for individ- 
ual isotopic chains. It was shown in (Jgl) that larger 



values of x always shift GTGR to higher excitation en- 
ergy but the shift caused by a small change of x does 
not result in the same energetic shift for all isotopes. 
On the other hand, the particle-particle-interaction (k) 
is known to enhance the ground-state correlations and 
leads to a redistribution of the calculated (3 strength 
(and is therefore related to the half- lives of the nuclei), 
which is commonly shifted to lower excitation energies 
([6|). Again, rather than going for a parameterization of 
K, the local best value for an isotopic chain that best 
reproduces the measured half-lives, was chosen in this 
work. Fixing the interaction constants, x f^nd k, in this 
way led to considerable improvement in the accuracy 
of the pn-QRPA calculation (0; @; S 0). For details of 
choosing the locally best values of the strength parame- 
ters see Ref. dsl). The chosen values of X were 0.27 MeV 
and 0.36 MeV for Kr and Sr isotopes, respectively. The 
corresponding k values were 0.045 Mev and 0.044 MeV, 
respectively. 

The deformation parameter was recently argued as 
an important parameter for QRPA calculations at par 
with the pairing parameter by Stetcu and Johnson (J18l ). 
To make the problem more challenging the nuclei in the 
mass region considered for this work are known to have 
coexistence of prolate and oblate shapes. For ^^Kr the 
deformation parameter was calculated using the mass 
formula of Ref. (|l9l ). The chosen value of -0.2647 is 



also in line with the findings of authors in Ref. (|l5l : l20f ) 
which favor an oblate shape for ^^Kr. For the case of 
^^Kr and ^^Sr. the experimentally adopted values of 
the deformation parameters, extracted by relating the 
measured energy of the first 2+ excited state with the 
quadrupole deformation, were available and taken from 
Raman et al. (|2lh . The deformation parameter for the 
case of ^^Sr was again taken from the mass formula of 
Ref. (Il9f ) in close agreement with the findings of Nacher 
and collaborators (|22l ) . Tabic [1] shows the chosen val- 
ues of nuclear deformations used in the model and the 
calculated half-lives for the waiting point nuclei. All 
half-lives are given in units of sec. Measured values 
are given in third column whereas the pn-QRPA cal- 
culated terrestrial half-lives of Kr and Sr isotopes are 
given in fourth column. The calculated terrestrial half- 
lives have dominant contribution from positron (Z?"*") 
decay rates. Since the chosen values of x and k repro- 
duced the measured half-lives fairly well, no quenching 
factor was introduced in the present calculation. The 
last column compares the pn-QRPA calculated stellar 
half-lives with the terrestrial ones under limiting con- 
ditions. At the lowest temperature considered in this 
work (T = 10^ K), excited parent states are not appre- 
ciably populated, while at the lowest density considered 
in this work (pYe = 10°'^ gcm~^), the continuum elec- 
tron density is quite low. It is to be noted that p stands 



for the baryon density whereas Ye is the ratio of the lep- 
ton number to the baryon number and their product 
pYe implies the electron density and is simply referred 
to as density throughout this manuscript. Except for 
electron capture, the stellar rates should be close to the 
terrestrial values. The exception arises because only 
continuum electron capture and no bound-states cap- 
ture is calculated. Again, for the stated physical condi- 
tions, the continuum electron capture rates are orders 
of magnitude smaller than the corresponding positron 
{I3~^) decay rates. As can be seen the comparison is 
satisfactory. Nuclear masses and Q-values required for 
the calculation were taken from the mass compilation 
of Audi et al. (|23[ ). The same set of model parameters 
were used to calculate electron capture and positron de- 
cay rates for these waiting point nuclei under rp-process 
conditions. 

The electron capture (ec) and positron decay (pd) 
rates of a transition from the z*'' state of the parent to 
the j*'' state of the daughter nucleus are given by 
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The value of D was taken to be 6295s (1241). B(F) and 
B(GT) are reduced transition probabilities of the Fermi 
and Gamow- Teller (GT) transitions, respectively. 
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Here (t*^ is the spin operator and i^ stands for the 
isospin raising and lowering operator with (gA/gv) = 
-1.254 (|25h. Details of the calculation of reduced tran- 
sition probabilities (for both ground state and excited 
states GT transitions) can be found in Ref. (J8|). 
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are the phase space integrals and are 



functions of stellar temperature (T), electron density 
(pYe) and Fermi energy (Ef) of the electrons. They are 
explicitly given by 






and by 
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Ill Eqs. (|4]) and (O, w is the total energy of the 
electron including its rest mass. Wm is the total /3-decay 
energy, 



nip — rrid + Ei — E. 
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where nip and Ei are masses and excitation energies of 
the parent nucleus, and m^ and Ej of the daughter nu- 
cleus, respectively. F(± Z,w) are the Fermi functions 
and were calculated according to the procedure adopted 
by Gove and Martin (J26[ ). G± are the Fermi-Dirac dis- 
tribution functions for positrons (electrons). 
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here E is the kinetic energy of the electrons and k is 
the Boltzmann constant. 

The total capture/decay rate per unit time per nu- 
cleus is finally given by 
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where Pi is the probability of occupation of parent ex- 
cited states and follows the normal Boltzmann distri- 
bution. After the calculation of all partial rates for the 
transition i — !■ j the summation was carried out over 
300 initial and final states and satisfactory convergence 
was achieved in the rate calculation. These 300 excited 
states were chosen with an appropriate bandwidth cov- 
ering energies up to 10 MeV in parent and 20 MeV in 
daughter nuclei. 

In order to ensure a reliable set of calculated rates 
three important steps were undertaken in this calcu- 
lation: (i) all available experimental (XUNDL) data 
were incorporated in the calculation; (ii) all available 
measured GT strength distribution were inserted in the 
calculation wherever the theory was missing them, and 
(iii) all excited state GT strength functions were micro- 
scopically calculated within the pn-QRPA framework. 
The theoretical excitation energy of a low-lying level 
was replaced by the observed excitation energy if this 
experimental level had a definite spin-parity assignment 
and was located within 0.5 MeV of the theoretically cal- 
culated level. If more than one theoretically calculated 
energy eigenvalue falls within 0.5 MeV of a measured 
level then the replacement was done for the level lying 
closest to the measured value. If there appeared a level 
without definite spin and/or parity assignment in the 
sequence of experimental levels from the ground state 
to higher excited levels, no replacement of theoretical 



levels with measured levels was done beyond this exci- 
tation energy. The calculated B(GT) value was then re- 
placed by the observed B(GT) value (if available) when 
the calculated energies of the parent and daughter lev- 
els were replaced by the experimental energies. Miss- 
ing measured states were inserted and inverse transi- 
tions (along with their log/t values) were also incorpo- 
rated in the calculation. If the relevant transition was 
not measured but the inverse transition was measured, 
the deduced observed B(GT) value using this inverse 
transition was used in this calculation. The measured 



B(GT+) strength distribution of "^^Kr from Ref. (jlj), 
of ^^Kr from Ref. Q and of ^®Sr from Ref. ^ were 
inserted in the calculation. As mentioned earlier the pn- 
QRPA model has the additional advantage that all ex- 
cited state GT strength functions are calculated in a mi- 
croscopic fashion and the so-called Brink's hypothesis 
is not assumed in this calculation. Brink's hypothesis 
states that GT strength distribution on excited states is 
identical to that from ground state, shifted only by the 
excitation energy of the state. Recent pn-QRPA calcu- 
lations (e.g. (|l2t Il3l )) have shown Brink's hypothesis 
to be a poor approximation to be used in calculation 
of stellar weak rates specially at high temperatures and 
densities. These steps ensured a reliable set of calcula- 
tion of ground and excited state GT strength distribu- 
tions to be later used for the calculation of weak rates 
for waiting point nuclei. 

Fig. [H shows the calculated BGT+ distributions for 
the Kr and Sr isotopes up to 20 MeV in daughter nuclei. 
As can be seen from Fig. [H the GT distributions are 
well fragmented and extend up to high excitation ener- 
gies in daughter nuclei. The model independent Ikeda 
sum rule was satisfied in the calculation. The discus- 
sion on the calculated weak rates follows in the next 
section. 



3 Thermal electron capture and positron decay 
rates of waiting point nuclei 

The calculated Fermi and GT strength distributions, 
along with the corresponding phase space calculations 
were inserted into Eq. ([T]) for the calculation of partial 
rates. The calculated partial rates were finally added 
through the recipe of Eq. ^ to get the total stellar 
electron capture and positron decay rates. The calcu- 
lated rates for Kr and Sr isotopes are shown in Fig. [2] 
and Fig. [3l respectively. 

The left panels of Fig. [2] and Fig. |3] depict the calcu- 
lated electron capture and positron decay rates, sepa- 
rately, as a function of stellar temperature and density. 
The electron capture rates are calculated at densities 



pYe = 10^ 10^ 10'^■^ 10'^ ana 10^ gcm-3. Tg gives the 
stellar temperature in units of 10^ K. The calculation 
is performed up to Tg = 30. It is to be noted that 
such high temperatures are far away from the relevant 
range of temperatures for the rp-process. Nevertheless 
the calculation was performed till Tg = 30 in order to 
study the evolution of rates with increasing stellar tem- 
peratures. The right columns in both figures show the 
sum of calculated electron capture and positron decay 
rates at the respective temperature and density points. 
The left panels show that the calculated positron de- 
cay rates are almost independent of stellar densities. 
It can be seen from the upper left panel of Fig. [2] 
that at low stellar temperatures (Tg < 4K) and den- 
sities (pYe < 10^'^ gcm""^), the positron decay rates of 
^^Kr are bigger than the electron capture rates by as 
much as a factor of 23. According to studies by au- 
thors in Ref. (jsl; y) the peak conditions for rp-process 
are in the vicinity of Tg = 1-3 K and pY^ = 10^ - 10^ 
gcm~^. It can be seen from the current calculation 
that electron capture and positron decay rates of ^^Kr 
contribute roughly equally to the total weak rates at 
rp-process conditions. However at high temperatures 
and densities the electron capture rates shoot up and 
are bigger by more than two orders of magnitude than 
the corresponding positron decay rates. For the case of 
^"^Kr the electron capture rates are always bigger than 
the positron decay rates (lower left panel of Fig. [2]). 
The electron capture rates are more than 1-2 orders 
of magnitude bigger than positron decay rates for rp- 
process conditions and the total rates are commanded 
by the electron capture rates (lower right panel). For 
still higher temperatures and densities the positron de- 
cay rates are almost 3 orders of magnitude smaller and 
can safely be neglected as compared to electron capture 
rates. 

Fig. [3] shows the calculated weak-interaction medi- 
ated rates for ^^-^^Sr. The upper panels are for ^^Sr and 
the lower for ^^Sr. Akin to the case of ^^Kr, at low stel- 
lar temperatures (Tg < 4.5K) and densities {pYe < 10^ 
gcni^"^), the positron decay rates dominate and are as 
much as a factor 17 bigger than the corresponding elec- 
tron capture rates. At high temperatures and densities 
the electron capture rates shoot exponentially due to 
increase in phase space factors and are more than two 
orders of magnitude bigger. For the rp-process condi- 
tions the situation is interesting. At density pY^ = 10^ 
gcm"-^ the positron decay rates are a factor 2-3 bigger, 
at pYe = 10^'^ gcm~'^ the two rates are of comparable 
magnitude while at pYg = 10^ gcm"-^ the electron cap- 
ture rates are a factor 4 bigger. At high temperatures 
and densities the total weak rates are driven by electron 
capture rates (upper right panel of Fig. [3]) . For the case 



of ^®Sr the positron decay rates are bigger than the elec- 
tron capture rates by as much as a factor of 7 at low 
temperatures (Tg < 3.5K) and density {pYe < 10^'^ 
gcm^'^) conditions. For the rp-process conditions the 
two rates are of equal magnitude at pY^ = 10^ gcm^'^ 
whereas the electron capture rates are bigger by a fac- 
tor of 3-10 at pYe = 10^-^ - 10'^ gcm-3. Once again at 
high temperatures the electron capture rates are bigger 
by more than two orders of magnitude as expected. 

The contribution of parent excited states to the to- 
tal electron capture and positron decay rates is shown 
in Table [2] The analysis was done at the selected stel- 
lar density of pY^ = 10^'^ gcm~^ (a typical value for 
rp-conditions) . Shown also in Table [2] are the ground 
state electron capture and positron decay rates and the 
ratios of the ground state capture rate to total elec- 
tron capture rate [Rec{G/T)) and ground state decay 
rate to total positron decay rate {Rpd{G/T)) for the 
waiting point nuclei. For the selected density scale no 
significant contribution comes to the total decay rates 
from parent excited states whereas at Tg = 30K the ex- 
cited states contribute roughly 95% to the total capture 
rates. It is to be noted that this analysis can change 
appreciably for higher stellar densities. 

Fig. [5] displays the half-lives (including both the 
positron decay and electron capture contributions), as a 
function of stellar temperatures at a selected density of 
pYe = 10^"'' gcm"'^ for Sr and Kr waiting point nuclei. 
Again the calculation is performed up to Tg — 30K for 
reasons mentioned above. The calculated half-lives de- 
crease as stellar temperature increases as expected since 
the electron capture rates increases substantially with 
increasing temperatures (see Fig. [2] and Fig. |3|). It is 
to be noted that for rp-process conditions the half-lives 
decrease very mildly. For higher temperature the fall 
in half-lives is rather steep. 

The calculated rates were also compared with the 
calculation of weak-interaction rates for Kr and Sr iso- 
topes by Sarriguren (J4|). The same author later ex- 
panded his calculation for waiting point nuclei extend- 
ing from Ni to Sn (J28h . Sarriguren used quasiparti- 
cle basis corresponding to a deformed self-consistent 
Skyrmc Hartec-Fock calculation with SLy4 (|2^) force. 
Pairing correlations were treated in BCS formalism. 
The spin-isospin interactions were considered both in 
particlc-holc and particle-particle channel and wave- 
functions were calculated using the Bohr-Mottelson fac- 
torization (|30l). Sarriguren used a quenching factor of 
0.55 in his calculation. The pn-QRPA calculated elec- 
tron capture rates are in good agreement with Sarrig- 
uren's rate except for the case of "^^Kr. The pn-QRPA 
positron decay rates are bigger specially at high tem- 
peratures by as much as a factor of 40. 



Tabic [3] compares the weak- interaction rates and cal- 
culated half-lives for Kr isotopes at pYe = 10^ gcm"'^. 
In the table Rec{N/S), Rpd{N/S) and Rhi{N/S) denote 
the ratio of pn-QRPA calculated electron capture (ec), 
positron decay (pd) and half-lives (hi), respectively, to 
those calculated by Sarriguren. It can be seen from 
Table |3] that the calculated electron capture rates on 
^^Kr are in reasonable agreement with those of Sar- 
riguren. Positron decay rates of ^^Kr are also in rea- 
sonable agreement except at Tg = lOK where reported 
rates are around a factor 5 bigger. Sarriguren did not 
calculate weak-interaction rates for Tg > lOK. The last 
column shows the ratio of the pn-QRPA calculated to- 
tal half-lives (including electron capture and positron 
decay) to those calculated by Sarriguren. The compar- 
ison of total half-lives for ^^Kr is correspondingly fairly 
good. Only at high temperatures is the pn-QRPA cal- 
culated half life roughly half that reported by Sarrig- 
uren. The origin of this difference might be connected 
to the fact that energy levels beyond 1 MeV were not 
considered in Sarriguren's calculation. These high-lying 
energy states have a finite occupation probability at 
high temperatures and resulted in a much bigger total 
weak rate. Much bigger differences are seen for the case 
of weak rates for ^^Kr. The reported electron capture 
rates on ^^Kr are a factor 3-7 bigger whereas positron 
decay rates are almost double at low temperatures and 
around a factor 30 bigger at Tg = lOK. The total weak 
rates are around an order of magnitude bigger and con- 
sequently the calculated half-lives are roughly an order 
of magnitude smaller. The reason for this difference 
could be attributed to the different calculation of the 
GT distributions in the two models. The pn-QRPA 
model calculated a total B(GT) value of 10.73, up to 20 
MeV, with the centroid lying around 5.10 MeV in ^^Br 
(see Fig. [1]). The pn-QRPA calculated total B(GT) 
strength value and/or placement of GT centroid can 
lead to enhancements in reported rates. 

A better comparison of reported calculation with 
that of Sarriguren is seen for the case of Sr isotopes (Ta- 
ble E]). For ^^Sr the electron capture and positron decay 
rates are in reasonable agreement. The calculated half- 
lives also agree fairly well except at Tg = lOK when the 
reported positron decay rate on ''^Sr is around a factor 
4 bigger. The comparison of electron capture rates on 
^*Sr is reasonable. The pn-QRPA reported positron de- 
cay rates are double the Sarriguren's rate at low tem- 
peratures and around factor 10 -40 bigger at higher 
temperatures. For the rp-process conditions the pn- 
QRPA calculated electron capture rates on ^^Sr (^*Sr) 
are slightly bigger (smaller) and positron decay rates for 
76gj. (78gj,-j g^j.g giig];i^;iy gmalkr (bigger) than the corre- 
sponding Sarriguren rates. The calculated half-lives of 



Sr isotopes are in reasonable agreement except at high 
temperatures when the half-life calculated by Sarrig- 
uren is around a factor 1.5 bigger for reasons already 
mentioned above. 



4 Conclusions 

To summarize weak-interaction rates for neutron- 
deficient Kr and Sr isotopes were calculated within the 
framework of pn-QRPA model. Proper choice of defor- 
mation parameter and smart choice of particle-particle 
and particle-hole interaction terms led to correct repro- 
duction of the measured half-lives. The same model 
parameters were used to calculate stellar weak rates for 
rp-process conditions. Whereas electron capture and 
positron decay rates, of other krypton and strontium 
isotopes in stellar matter, can be calculated using the 
same model parameters as stated here, a case by case 
fitting procedure is required for other isotopic chains 
within the framework of pn-QRPA model. Measured 
data were inserted wherever available to increase the 
reliability of calculated rates. The pn-QRPA reported 
electron capture rates on ^^Kr are around a factor 6 
bigger than those calculated by Sarriguren under rp- 
conditions. The present calculation clearly shows that 
the electron capture rates are at least of similar magni- 
tude as the competing positron decay rates under rp- 
process conditions. The calculation provide concrete 
evidence that for typical rp-process conditions, Tg = 
1-3 K and pYp = 10^ — 10^ gcm^-^, electron capture 
rates on ^^Kr and ^^Sr are of same order of magnitude 
as the corresponding positron decay rates. On the other 
hand electron capture rates on ^^Sr and ^^Kr are 1-2 
orders of magnitude bigger than the respective positron 
decay rates. The study reconfirms the conclusion made 
by Sarriguren that electron capture rates form an in- 
tegral part of weak-interaction mediated rates under 
rp-process conditions and should not be neglected in 
nuclear reaction network calculations as done in past 
(e.g. (0)). 
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Table 1 Comparison of measured and pn-QRPA calculated half-lives in units of sec. Nuclear deformations are given in 
second column. The third and fourth columns give the measured and calculated terrestrial half-lives, respectively. The fifth 
column gives the limiting values of calculated stellar half-lives at density lO'^'^gcm'^ and temperature 10^ K. 



Nucleus 


Deformation 


Ti/2(exp) 


ri/2(terr) 


Ti 


/2(0.5,0.01) 


72Kr 


-0.2647 


17.2 


17.2 




16.7 


74Kr 


-H0.3870 


690 


683.8 




665.7 


™Sr 


-H0.4074 


8.9 


8.8 




8.53 


78Sr 


-H0.4340 


150 


154.7 




142.95 



Table 2 The ground state electron and positron decay rates, Xec{G), Xpd{G), respectively, for ''^''^'^Kr and ''^'^'''^'Sr in 
units of sec~^. Given also are the ratios of the ground state capture and decay rates to total rate, Rec{G/T), Rpd{G/T), 
respectively. The first column gives the corresponding values of stellar density, pYe {gcm~^), and temperature, Tg (in units 
of 10^ K), respectively. 







72Kr 






(pYe,T9) 


Aec(G) 


Rec(G/T) 


Apd(G) 


Rpd(G/T) 


(lo6■^o.ol) 


3.57E-02 


l.OOE+00 


4.15E-02 


l.OOE+00 


(106-5,1) 


3.59E-02 


l.OOE+00 


4.15E-02 


l.OOE+00 


(106-5,1.5) 


3.62E-02 


l.OOE+00 


4.15E-02 


l.OOE+00 


(10^-5,2) 


3.67E-02 


l.OOE+00 


4.15E-02 


l.OOE+00 


(10^-5,2.5) 


3.75E-02 


l.OOE+00 


4.15E-02 


l.OOE+00 


(106-5,3) 


3.91E-02 


l.OOE+00 


4.15E-02 


l.OOE+00 


(106-5,30) 


1.48E+01 


2.83E-02 


2.24E-03 


6.69E-04 


(106-5,0.01) 


4.15E-02 


74^1. 

l.OOE+00 


1.04E-03 


l.OOE+00 


(106-5,1) 


4.25E-02 


l.OOE+00 


1.04E-03 


l.OOE+00 


(106-5,1.5) 


4.36E-02 


l.OOE+00 


1.04E-03 


l.OOE+00 


(106-5,2) 


4.52E-02 


l.OOE+00 


1.04E-03 


l.OOE+00 


(106-5,2.5) 


4.74E-02 


l.OOE+00 


1.04E-03 


l.OOE+00 


(106-5,3) 


5.05E-02 


l.OOE+00 


1.04E-03 


9.98E-01 


(106-5,30) 


1.16E+01 


5.37E-02 


3.26E-05 


9.61E-05 


(106-5,0.01) 


9.61E-02 


l.OOE+00 


8.12E-02 


l.OOE+00 


(106-5,1) 


9.69E-02 


l.OOE+00 


8.12E-02 


l.OOE+00 


(106-5,1.5) 


9.79E-02 


l.OOE+00 


8.12E-02 


l.OOE+00 


(106-5,2) 


9.96E-02 


l.OOE+00 


8.12E-02 


l.OOE+00 


(106-5,2.5) 


1.02E-01 


l.OOE+00 


8.12E-02 


l.OOE+00 


(106-5,3) 


1.07E-01 


l.OOE+00 


8.12E-02 


l.OOE+00 


(106-5,30) 


2.36E+01 


3.93E-02 


4.42E-03 


1.23E-03 


(106-5,0.01) 


1.42E-02 


l.OOE+00 


4.85E-03 


l.OOE+00 


(106-5,1) 


1.43E-02 


l.OOE+00 


4.85E-03 


l.OOE+00 


(106-5,1.5) 


1.45E-02 


l.OOE+00 


4.85E-03 


l.OOE+00 


(106-5,2) 


1.48E-02 


l.OOE+00 


4.85E-03 


l.OOE+00 


(106-5,2.5) 


1.53E-02 


l.OOE+00 


4.85E-03 


l.OOE+00 


(106-5,3) 


1.60E-02 


l.OOE+00 


4.84E-03 


9.98E-01 


(106-5,30) 


3.59E+00 


2.60E-02 


1.41E-04 


1.77E-04 



Table 3 Comparison of electron capture ( 
Kr isotopes. Number in fourth and seventh 
The last column shows the ratio of the pn 
are compared at stellar density pYe ~ 10^ 



ec) and positron 
column gives the 
QRPA calculated 



decay (pd) rates with those calculated by Sarriguren (|j) for 
ratio of reported rates to the rates calculated by Sarriguren. 
total half-lives to those calculated by Sarriguren. The rates 



gem 







ec 


72 


Kr 


pd 






Tg 


Nabi 


Sarriguren 


B:e.{N/S) 


Nabi 


Sarriguren 


Rpd{N/S) 


B,u{N/S) 


0.01 


1.29E-02 


1.62E-02 


7.97E-01 


4.15E-02 


3.47E-02 


1.20E+00 


9.35E-01 


1 


1.23E-02 


1.54E-02 


7.99E-01 


4.15E-02 


3.48E-02 


1.19E+00 


9.33E-01 


2 


1.19E-02 


1.55E-02 


7.66E-01 


4.15E-02 


3.99E-02 


1.04E+00 


1.04E+00 


3 


1.37E-02 


1.92E-02 


7.12E-01 


4.15E-02 


5.07E-02 


8.19E-01 


1.27E+00 


10 


1.33E+00 


9.61E-01 


1.38E+00 


3.38E-01 


7.52E-02 


4.49E+00 


6.21E-01 


0.01 


1.43E-02 


2.20E-03 


6.49E+00 


Kr 

1.04E-03 


5.48E-04 


1.90E+00 


1.79E-01 


1 


1.40E-02 


2.20E-03 


6.36E+00 


1.04E-03 


5.79E-04 


1.80E+00 


1.85E-01 


2 


1.43E-02 


2.66E-03 


5.38E-I-00 


1.04E-03 


8.78E-04 


1.18E+00 


2.30E-01 


3 


1.74E-02 


3.65E-03 


4.76E+00 


1.03E-03 


1.12E-03 


9.20E-01 


2.59E-01 


10 


1.04E+00 


3.32E-01 


3.13E+00 


4.44E-02 


1.44E-03 


3.07E+01 


3.08E-01 



Table 4 Same as Table [3] but for Sr isotopes. 







ec 


76 


Sr 


pd 






Tg 


Nabi 


Sarriguren 


Reo{N/S) 


Nabi 


Sarriguren 


Rpd{N/S) 


Rhi{N/S) 


0.01 


3.45E-02 


3.04E-02 


1.13E+00 


8.12E-02 


9.28E-02 


8.7.5E-01 


1.06E+00 


1 


3.31E-02 


2.99E-02 


l.llE+00 


8.12E-02 


9.95E-02 


8.16E-01 


1.13E+00 


2 


3.22E-02 


3.04E-02 


1.06E+00 


8.12E-02 


1.14E-01 


7.10E-01 


1.28E+00 


3 


3.72E-02 


3.55E-02 


1.05E+00 


8.11E-02 


1.22E-01 


6.63E-01 


1.33E+00 


10 


2.20E+00 


1.55E+00 


1.42E+00 


4.73E-01 


1.37E-01 


3.46E+00 


6.33E-01 


0.01 


5.09E-03 


5.63E-03 


9.04E-01 


Sr 

4.85E-03 


2.46E-03 


1.97E+00 


8.14E-01 


1 


4.88E-03 


5.66E-03 


8.63E-01 


4.85E-03 


2.75E-03 


1.76E+00 


8.64E-01 


2 


4.78E-03 


6.07E-03 


7.87E-01 


4.85E-03 


3.35E-03 


1.45E+00 


9.78E-01 


3 


5.59E-03 


7.48E-03 


7.47E-01 


4.85E-03 


3.59E-03 


1.35E+00 


1.06E+00 


10 


5.72E-01 


5.40E-01 


1.06E+00 


1.79E-01 


4.60E-03 


3.90E+01 


7.25E-01 
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Fig. 1 Cumulated Gamow- Teller strength distributions B(GT+) for Kr and Sr isotopes. The abscissa represents energy 
in units of MeV in daughter nucleus. 
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Fig. 2 (Color online) Weak-interaction mediated rates for Kr isotopes as a function of stellar temperature and density. 
The left panels show the electron capture and positron decay rates whereas the right panels show the combined rates. All 
rates are given in units of s~^. Densities are given in units of gcm~^ and temperatures in units of 10*' K. 
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Fig. 3 (Color online) Same as Fig. [2] but for Sr isotopes. 
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Fig. 4 (Color online) Half-lives for Kr and Sr isotopes as a function of stellar temperature calculated at a density of 



